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SUMMARY 

A  staple  constitutive  law  Is  proposed  for  the  description  of  a  ceramic 
composite  which  undergoes  stress  Induced  martensitic  transformation*  This  law 
is  used  in  finite  element  calculations  to  investigate  the  shear  effect  on  the 
transformation  zone  near  a  crack  tip.  A  formula  describing  the  Btress 
intensity  factor  change  due  to  the  shear  contribution  of  the  transformation  is 
given.  Significant  loss  of  toughness  Is  observed  in  the  case  of  a  stationary 
crack  and  is  attributed  entirely  to  the  shear  component  of  the  transfor¬ 


mation. 


about  no 


As  the 


crack  grows,  the  wake  of  the  transformed  material  left  behind  the  crack 
constitutes  a  source  of  toughening.  This  toughening  is  due  to  both  dilatancy 
and  shear  in  the  phase  change  and  rises  to  a  maximum  level  just  after  a 
propagation  comparable  with  the  zone  height.  Finally,  it  is  shown  that  the 
shear  component  can  be  important  when  prediction  of  the  fracture  toughness  of 
the  transformation  toughened  ceramics  are  made. 
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1.  INTRODUCTION 


Transformation  toughening  is  one  of  the  mechanisms  available  to  overcome 
the  Inherent  brittleness  of  ceramics*  It  is  a  phenomenon  applicable  to 
ceramic  matrices  in  which  Zirconia  (Zr02>  and  perhaps  some  other  materials  can 
be  incorporated*  To  date  it  has  been  studied  quite  thoroughly  [1-13]. 

An  optimally-fabricated  [13]  partially  stabilized  zirconia  (FSZ)  is  a 
two-phase  ceramic*  Its  microstructure  [5]  at  room  temperature  consists  of:  _a 
cubic  matrix  which  is  a  high  solute  content  Zirconia  'alloy’  containing  one  of 
the  stabilizers  MgO  CaO,  Y2O3  or  any  of  the  rare  earth  oxides:  and  fine 
coherent  metastable  tetragonal  precipitates  of  low  solute  content  Zirconia 
phase  inside  the  cubic  matrix*  The  stress  Induced  martensitic  transformation 
of  those  meta8tably  retained  tetragonal  particles  to  monocllnlc  symmetry  in 
the  stress  field  of  a  crack  tip  is  the  mechanism  responsible  for  the  enhanced 
toughness  observed  experimentally  in  PSZ.  The  phenomenon  of  transformation 
toughening  is  also  observed  in  the  AI2O3  and  Zr(>2  system.  The  pure  tetrogonal 
Zr02  particles  are  retained  metastably  in  the  AI2O3  matrix  by  pressure. 

There  are  two  methods  of  analysis  with  regard  to  the  phenomenon  of 
toughening.  The  first  Incorporates  the  energy  changes  accompanying  the 
transformation*  The  second  concerns  the  stress  intensity  factor  (SIF)  changes 
that  take  place  during  transformation. 

Quantitative  analysis  from  a  continuum  mechanics  viewpoint,  based  only  on 
the  dllatational  component  of  the  transformation  with  regard  to  the  SIF 
reduction  was  done  first  by  McMeeking  and  Evans  [19]  and  Budiansky,  Hutchinson 
and  Lambropoulos  [20].  The  analysis  of  the  latter  workers  is  based  on  a 
constitutive  relation  between  the  mean  stress  and  the  dilatation  for  the 
composite  ceramic.  In  both  works  mentioned,  the  predicted  toughening  is 
comparable  with  the  experimental  data,  however,  it  underestimates  them. 
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Recently,  Lambropoulos  [21]  has  suggested  a  promising  constitutive  law 
for  the  composite  including  both  parts  of  the  transformation,  dilatant  and 
shear.  The  treatment  of  the  shear  effect  takes  into  consideration  particle 
size  and  orientation*  The  results  predicted  for  spherical  particles  are  quite 
well  in  agreement  with  experimental  observations  but  are  based  on  a 
transformation  zone  size  and  shape  estimated  from  the  standard  crack  tip 
singular  elastic  field.  That  is,  he  assumes  that  the  zone  shape  Is  the  same 
as  regions  in  the  unperturbed  elastic  solution  in  which  the  transformation 
criterion  is  met  or  exceeded.  Changes  in  stress  due  to  the  transformation  are 
not  addressed  in  detail. 

In  this  paper,  we  shall  study  the  shear  effect’s  influence  on  the 
enhanced  toughness  through  numerical  calculations  by  means  of  a  simple  stress- 
strain  relation  for  the  composite.  The  model  is  based  on  a  constitutive 
relation  along  the  lines  Introduced  by  Budiansky  et  al.  [20]  and  the  condition 
for  transformation  is  dominated  by  the  dilatatlonal  component.  When  the 
hydrostatic  stress  reaches  a  critical  level,  the  transformation  takes  place. 
This  model  is  not  entirely  satisfactory  as  the  shear  effect  is  bound  to 
Influence  the  critical  state  for  transformation.  However,  we  regard  this 
paper  with  the  critical  state  determined  solely  by  the  hydrostatic  stress  as  a 
first  step  towards  a  complete  theory.  The  analysis  is  carried  out  around  the 
crack  tip  of  a  long  crack  which  is  Imbedded  in  a  composite  material  rich  in 
tetragonal  particles  whose  presence  in  the  matrix  is  defined  by  a  volume 
concentration  Vf.  As  such,  the  composite  can  be  modelled  as  a  continuum  of 
transformable  material. 

In  section  2  we  discuss  some  aspects  of  the  shear  part  of  the 
transformation  which  are  helpful  in  comprehending  the  nature  of  the  shear 
strain  that  the  transforming  particles  undergo.  In  section  3  we  propose  a 


constltudve  law  for  the  composite*  In  section  4  we  derive  a  formula  for  the 
SIF  change  due  to  the  shear  component  of  the  transformation.  Based  on  that 
formula,  we  make  first  estimates  of  the  shearing  contribution  to  SIF  change  by 
using  the  transformation  zone  derived  from  the  unperturbed  elastic  solution 
[19].  In  section  5  we  formulate  the  boundary  value  problem  for  the  stationary 
crack  and  solve  it  by  means  of  the  finite  element  method  in  section  6.  As  a 
result  the  transformation  zone  shape  and  size  are  estimated.  In  section  7  the 
estimated  zone  and  the  fracture  toughness  calculations  for  the  stationary  and 
for  the  propagating  crack  are  presented.  As  it  has  already  been  proven  In  the 
past,  the  toughness  Increase  is  due  to  the  transformed  particles  left  in  the 
wake  of  the  crack  tip.  In  sections  8  and  9  the  discussion  associated  with  the 
model  results  and  the  closure  are  presented  respectively. 

2.  STRESS  INDUCED  TRANSFORMATION 

The  aim  of  this  section  is  to  introduce  basic  features  of  the  stress 
Induced  transformation  so  that  the  shear  component  can  be  understood.  The 
transformation  is  martensitic  and  has  been  discussed  extensively  elsewhere 
[1,16,17,18].  It  Involves  a  change  from  tetragonal  to  monocllnlc  symmetry  in 
particles  in  the  composite  ceramic.  In  the  situation  of  interest  to  us,  the 
transformation  is  Induced  by  critical  conditions  of  stress.  A  component  of 
the  transformation  is  a  dilatation  and  if  the  particles  were  unconstrained 
there  would  also  be  a  substantial  shear  contribution.  If  we  assume  that  the 
process  is  driven  by  the  reduction  of  free  energy  [1,7,13]  then  we  deduce  that 
the  shear  strain  component  of  the  transformation  would  align  itself  to 
maximize  the  work  done  by  the  loads  applied  to  the  particle  [7,13].  In 
addition,  the  critical  state  for  transformation  would  arise  just  when  the 
applied  loads  are  capable  of  delivering  sufficient  energy  to  the  system  to 
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compensate  for  the  increase  In  the  internal  energy  [7,13).  This  implies  a 
transformation  criterion  involving  some  combination  of  hydrostatic  and 
devlatorlc  stress.  This  issue  has  been  addressed  by  Lambropoulos  [21]. 

It  has  been  observed  that  the  situation  is  more  complicated  when  the 
transforming  particles  are  constrained  in  the  composite  matrix  [12].  The 
particles  are  capable  of  twinning  or  undergoing  some  similar  mode  of 
deformation  during  transformation.  The  twins  form  in  such  a  way  that  the 
average  shear  strain  after  transformation  can  be  quite  small  compared  with  the 
potential  unconstrained  shear  transformation  strain.  In  addition,  the  final 
strain  in  the  transformed  particle  differs  because  of  the  constraint  of  the 
surrounding  matrix.  There  is  not  yet  a  comprehensive  theory  that  accounts  for 
transformation  and  twinning.  However,  it  can  be  hypothesized  that  the 
orientation  of  the  net  shear  strain  that  results  from  transformation  and 
twinning  will  be  aligned  with  the  maximum  shear  stress  applied  by  the  matrix 
to  the  particle.  This  will  tend  to  maximize  the  external  energy  absorption 
during  transformation  and  suggests  a  critical  state  for  the  transformation 
based  on  strain  energy.  Lambropoulos  [21]  has  developed  a  constitutive  law 
for  the  constrained  transformation  along  the  above  lines  combining  the  effects 
of  hydrostatic  and  devlatorlc  stress  and  accounting  for  dilatational  and  shear 
transformation. 

In  this  paper  we  shall  use  a  simpler  law  as  a  first  step  towards  studying 
the  interactions  of  shear  transformation  with  the  crack  tip.  The  process  will 
be  considered  to  take  place  on  a  continuum  scale  and  the  description  of  the 
constitutive  law  applies  to  the  composite  ceramic.  The  criterion  for 
transformation  will  be  taken  to  be  the  achievement  of  a  critical  average 
hydrostatic  stress  in  the  composite.  This  neglects  the  contribution  made  by 
the  shear  component  to  the  work  absorbed  during  transformation.  However,  the 


residual  shear  strains  In  constrained  transformation  particles  are  known  to  be 


small  and  the  shear  strain  in  an  unconstrained  composite  element  would  be 
correspondingly  small.  The  transformation  trill  involve  a  deviatorlc  component 
as  well  as  a  dilatatlonal  contribution.  The  deviatorlc  part  of  the 
transformation  strain  for  an  unconstrained  element  of  the  composite  ceramic 
will  be  taken  to  be  proportional  to  the  dilatatlonal  component  of  the 
unconstrained  transformation.  All  calculations  presented  in  this  paper  will 
be  for  plane  strain  situations.  In  that  case  there  is  a  fixed  ratio  between 
the  shear  strain  and  the  volumetric  strain  in  the  transformation.  The 
orientation  of  the  shear  strain  will  be  taken  as  that  of  the  maximum  shear 
stress  when  transformation  commences  (fig.  1).  Once  the  transformation  has 
taken  place  this  orientation  will  be  locked  in  so  that  changes  of  direction  of 
the  maximum  shear  stress  will  not  cause  rotation  of  the  shear  contribution  of 
the  transformation. 

Finally,  the  transformation  will  be  supercritical  in  the  terminology  of 
Budiansky  et  al.  [20].  That  is,  at  critical  state,  the  material  transforms 
completely  without  the  existence  of  a  partially  transformed  state.  The 
details  of  the  constitutive  law  are  given  in  the  next  section. 


3.  THE  CONSTITUTIVE  LAW 


The  transformation  of  the  composite  occurs  due  to  the  martensitic 
transformation  of  the  particles.  It  takes  place  when  the  macroscopic  average 
stress  in  the  composite  is  such  that 


where  o 

m 


a 

m 


is  the  mean  stress, 


is  a  critical  value  and 


(1) 

°ij  is 


the  macroscopic  average  stress  tensor  at  a  point  in  the  composite 


The  unconstrained  transformation  strain  of  the  composite  is  t i j 


which 


is  partly  dilatational  and  partly  deviatoric  in  general*  It  can  be  estimated 

to  be  the  volume  fraction  Vf  of  particles  in  the  composite  times  the 

unconstrained  transformation  strain  of  the  transforming  particles  less  the 

deviatoric  strain  nullified  by  twinning  [20].  The  supercriticality  of  the 

transformation  implies  that  there  is  no  transformation  as  long  as  0m  <  oJjjr 

and  complete  transformation  occurs  if  a  >  a  .  The  transformation  will  be 
r  mm 

assumed  to  be  effectively  irreversible  in  the  conditions  prevailing  as 
observed  in  experiment  [19] «  The  amount  of  the  deviatoric  strain  arising 
during  transformation  is  proportional  to  the  amount  of  dilatation.  We  will 
carry  out  calculations  for  a  variety  of  values  of  this  ratio.  The  principal 
axes  of  the  deviatoric  strain  will  be  taken  to  depend  on  the  state  of  stress 
at  transformation  as  discussed  below. 

The  stress  which  arises  [22]  in  a  constrained  element  of  the  composite  is 


°ij  "  cijk*  (ekA  "  ekP 


(2) 


where  e^j  is  the  final  strain  and  *8  t'le  ten8or  of  linear  elastic 

moduli  of  the  composite  material.  For  an  isotropic  composite 


°ij  "  (eij  ”  eij*  +  B(ekk  "  £kk)6ij 


(3) 


where  p  is  the  shear  modulus,  B  is  the  bulk  modulus,  e^j  *  tjj  -  y  6 


kk  ~ij 


is  the  deviatoric  strain,  djj  Is  the  Kronecker  delta  and  e1  i8  the 


dilatational  part  of  the  unconstrained  transformation  of  the  composite  (e 


ekk^ 

The  calculations  we  have  carried  out  are  for  plane  strain.  In  this  case 


ezz  ”  0  and  we  assume  that  ez z  ■  0  as  there  are  no  macroscopic  transverse 
shear  strains.  As  a  consequence 


v(o  +o  )  -  t  1*0  +  v)  e 
XX  yy  3 


(4) 


o 

zz 

where  v  Is  the  poisson's  ratio* 
written  as 


The  deviatoric  transformation  strain  can  be 


eT  -  -  eT  -  -  i-  sin2Q,  eT  -  i-  cos2Q 
xx  yy  2  xy  2 

where  Q  is  the  angle  between  the  x-axis  and  a  principal  axis  of  the 
transformation  shear.  Thus,  the  transformation  is  determined  by  3  parameters 
eT,  yT  and  Q.  We  will  assume  that  Q  is  coincident  with  the  angle  to  the 
principal  axes  of  shear  stress  in  the  macroscopic  composite  state  of  stress  at 
the  Instant  of  transformation  as  shown  in  fig.  1.  Calculations  are  carried 
out  for  a  variety  of  ratios  \  m  yT/eT. 

4.  ESTIMATES  OF  THE  SHEARING  EFFECT  AT  THE  CRACK  TIP 

Before  proceeding  to  somewhat  rigorous  numerical  calculations,  we  shall 
consider  some  approximate  results  for  shearing  transformation  at  the  crack 
tip.  The  shape  of  the  zone  of  transformed  material  at  the  crack  tip  is 
determined  by  the  Interaction  of  the  stresses  generated  by  the  applied  load 
and  those  generated  by  constraints  on  the  transformed  zone.  If  we  neglect  the 
latter,  we  can  estimate  the  zone  shape  as  the  locus  of  points  at  critical 
state  in  the  unperturbed  linear  elastic  solution  at  the  crack  tip.  This 
proves  to  be  quite  an  accurate  estimate  of  shape  for  the  case  of  small  scale 
dllatant  transformation  [20].  As  we  have  approximated  the  critical  state  as 
one  that  depends  only  on  hydrostatic  stress,  the  zone  shape  will  be  given  by 
the  locus  of  points  of  equal  hydrostatic  stress.  We  will  restrict  ourselves 
to  small  scale  transformation,  so  that  the  stresses  of  interest  are  given  by 
the  singular  elastic  stresses  at  the  crack  tip  due  to  the  applied  load 


VvV 


where  K^  is  the  stress  intensity  due  to  the  applied  loads  causing  tensile 
opening  (Mode  I)  of  the  tip,  (r,9)  are  polar  coordinates  measured  from  the 


crack  tip  and  f^j  is  a  given  function  which  can  be  found  for  example  in  the 
article  by  Rice  [24],  The  shape  that  results  for  a  stationary  crack  is  shown 
in  fig,  2. 

Consider  now,  the  material  which  transforms  inside  the  zone.  If  it  were 
not  constrained  by  the  material  outside  the  zone,  a  certain  change  of  shape 
would  result.  Tractions  can  be  applied  to  the  perimeter  of  this  region  of 
material  to  return  it  exactly  to  the  shape  of  the  zone  prior  to  transfor¬ 
mation.  After  the  material  is  inserted  into  the  crack  tip  location,  the 
tractions  can  be  removed  to  give  the  final  state.  However,  the  forces 
required  to  nullify  the  constraining  tractions  will  produce  a  change  in 

the  SIF  at  the  crack  tip.  It  Is  this  change  of  stress  intensity  AK  which  is 
of  interest. 

As  discussed  by  McMeeking  and  Evans  [19]  the  change  of  stress  intensity 
is  given  by 


AK  =*  J  Tc  •  h  ds  (6) 

ST 

where  _h_  is  the  weight  function  [25]  whose  form  is  stated  in  McMeeking  and 
Evans  [19]  and  used  with  the  assumption  that  tjj  is  homogeneous  in  A^>  ,  the 
transforming  area  which  has  perimeter  S-p.  We  shall  consider  now 
transformations  which  are  inhomogeneous  With  the  area  Aj  removed 

from  constraint  of  the  surrounding  area,  the  displacements  and  strains  [22] 

C  C 

that  result  will  be  uf  and  £  respectively  due  to  transformation  and  self- 
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constraint  within  Af.  This  will  give  rise  to  an  elastic  stress  such  that 


°h  ■  C1JU  <4  -  £u> 


"l  “lj  '  0  on  ST 


where  n  is  the  unit  outward  normal  to  ST  and  a .  .  .  m  t — —  with  x. 

—  1  ij.i  1 

denoting  position  in  a  fixed  cartesian  coordinate  system.  Tractions  are  now 

R  R 

applied  around  Aj  producing  further  displacements  uy  and  strains  £  so 
that  Aq>  has  its  original  shape  prior  to  transformation.  As  a  result 


R  .  /eR  ,  S  . 

ij  CijkX  (  kA  kJl  “  kA) 


o  =»  0 

ij.i 


R  .  S  n 
ui  +  ui  ’  0 


on  S  . 


Using  we  can  express  _TC  as  follows 


Thus,  equation  (6)  becomes 


Tj  “  ~  "i  °1J 
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As  noted  by  Rice  [25],  can  be  treated  as  a  displacement  and  so  by  the 


principle  of  virtual  work 


“  -  -  I  ^  hlfl  (9) 

T 

Substitution  of  (8a)  into  (9)  furnishes 

AK  - J  (ejj  -  ejj  -  tjj)  CljU  hl  k  dA  (10) 

X 

As  further  noted  by  Rice  [25]  ^ifk  *s  a  3tre88  field  in  equilibrium 

in  At  because  h^  is  related  to  the  difference  between  two  displacement 
fields  each  possessing  equilibrium  stress  fields.  This  means  that  the  virtual 
work  principle  can  be  used  to  show 

J  +  ejj)  CljkJL  hA  k  dAT  -  J  ni  hA>k(u*  +  u^)  dS 

»  0  due  to  (8c) 

This  means  that  equation  (10)  is  written  as 


AK  *  J  Eij  CijkJl  hA,k  ^  (11) 

For  isotropic  material,  equation  (11)  becomes 


AK 

-  AKd  +  AKS 

(12a) 

where 

4kd  ■*  J  *t  \,k  «* 

*T 

(dilatational  contribution) 

(12b) 

AKS  -  2d  J  .Jj  ht  j  dA 

(shearing  contribution) 

(12c) 

lv)  experimentally  measured  value 

AK  •  -  2.3  MPa  /m  (implies  X  =  3) 

The  preceding  toughness  enhancement  calculations  should  be  regarded 
tentative  because  of  the  uncertainties  involved  in  the  material  parameters. 
Nevertheless,  it  is  clear  that  the  predicted  toughness  increase  may  be  in 
agreement  with  the  experimental  data  when  the  shearing  contribution  is 
considered  in  the  transformation  mechanism  analysis.  An  important  role  is 
played  by  the  parameter  X  which  determines  the  amount  of  the  shearing 
contribution.  We  have  arbitrarily  used  1  because  we  have  no  information  on 
what  values  may  be  realistic.  Alternately,  we  can  compute  a  value  of  X 
necesssary  to  bring  our  estimates  into  agreement  with  the  experimental  values. 
These  are  the  implied  values  of  X  listed.  Some  of  these  values  are  rather 
large.  However,  the  amount  of  shear  strain  during  transformation  is  typically 
quite  large  compared  to  the  dilatation.  Even  if  substantial  amounts  of  this 
are  nullified  by  twinning,  this  could  still  leave  values  of  X  of  the  order 
Implied. 

9.  CLOSURE 

For  further  and  more  elaborate  treatment  of  the  shear  effect  one  may  need 
to  incorporate  size  and  particle  orientation  effects  via  a  more  realistic 
constitutive  law  for  the  composite.  This  must  include  a  consideration  of  the 
effect  of  twinning  during  the  transformation.  However,  our  simple 
calculations  have  indicated  that  shear  strain  effects  may  be  significant  and 
that  further  work  on  the  mechanics  of  the  phenomenon  may  be  profitable. 
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AK  -  -  1.70  MPa  /m 
iv)  experimentally  measured  value 

AK  -  -  2.4  MPa  /m  (Implies  X  -  1.7) 

T 

For  another  system  whose  the  parameters  are  {21]:  E  *  470  GPa,  e  *  0.04  w 
5.10“^ni,  Vf  ■  0.3  one  has 

i)  with  purely  dilatant  transformation 

AK  -  -  3.96  MPa  /m 

ii)  with  both  parts  of  the  transformation  and  X  *  1 

AK  *  -  6.30  MPa  v'm 

iii)  Lambropoulos'  result  (both  parts  of  the  transformation,  spherical 
particles  and  Kjq  -  6  MPa  /m) 

AK  -  -  4.5  MPa  /m 

iv)  experimentally  measured  value 

AK  -  -  6  MPa  /m  (implies  X  =  0.9) 

In  the  case  of  the  MgO  partially  stabilized  Zlrconla  whose  parameters  are 
[19]:  E  -  200  GPa,  eT  *  0.058,  w  -  6. 10~6m,  vf  *  0.3  one  has: 

i)  with  purely  dilatant  transformation 

AK  -  -0.85  MPa  /m 

ii)  with  both  parts  of  the  transformation  and  X  *  1 

AK  -  -  1.35  MPa  /m 

iii)  Lambropoulos '  result  (both  parts  of  the  tranformation,  spherical 
particles  and  Kj^  ■  1.97  MPa  /m) 


AK  -  -  1.04  MPa  /m 


sufficient  In  general  to  consider  only  the  effect  of  very  long  wakes  that 
produce  an  asymptotic  value  for  AK. 


Por  X  less  than  1  the  toughness  enhancement  is  not  as  large  as  that 
given  by  equation  (18).  In  table  5  we  present  the  greatest  toughness  enhance¬ 
ments  found  for  several  X  values  and  the  amount  of  crack  advance  for  which 
they  were  obtained.  This  table  shows  that  the  toughness  enhancement  gets 
larger  when  X  gets  larger.  It  should  be  mentioned  that  the  values  for  AK 
tend  to  certain  asymptotic  values  as  Aa/w  -*■  «  but  these  are  less  Important 
than  the  maxima  presented  in  table  5.  Finally,  we  should  mention  the  fact 
that  the  results  of  table  5  fit  the  curve  for  v  ■  0.25  of  figure  4.  This 
means  that  the  results  taken  in  section  4  and  those  after  having  solved  the 
boundary  value  problem  of  the  stationary  crack  are  the  same. 

Applications 

Next,  we  proceed  to  see  how  the  maximum  toughness  enhancement  results 
compare  with  the  experimental  data.  We  shall  consider  v  «  0.3  even  though 
our  shearing  contribution  results  have  been  taken  by  considering  v  »  0.25. 
Figure  4  indicates  that  no  significant  difference  results.  In  the  case  of  the 
AI2O3  toughened  Zlrconla  whose  material  and  transformation  parameters  are 
[21]:  E  *  315  GPa,  eT  ■  0.04,  w  -  10”^m,  v^  -  0.03  one  has 
i)  with  purely  dllatant  transformation 

AK  -  -  1.19  Mpa  /m 

ii)  with  both  parts  of  the  transformation  and  X  -  1 

AK  -  -  1.89  MPa  /m 

ill)  Lambropoulos '  result  (both  part  of  the  transformation,  spherical 
particles  and  KjC  -  5  MPa  /m) 


based  on  equation  (16).  Such  "R  curves’*  are  those  shown  In  figures  12  and  13 
for  X  =  0  and  X  *  1  respectively.  Observing  those  curves  and  the  tables  2 
through  4  we  deduce  that  the  wake  of  the  transformed  material  left  behind  the 


propagating  crack  is  a  source  of  fracture  toughness  enhancement.  First, 
figure  12  shows  that  the  toughness  enhancement  found  for  X  *  0  coincides 
with  that  found  by  McMeeking  and  Evans  [19].  This  coincidence  becomes  more 
pronounced  for  large  values  of  crack  advance. 

The  important  result  of  this  paper  is  the  AKg  component  that  behaves  as 
for  examaple  figure  13  shows  for  X  =  1.  As  mentioned  previously,  the  crack 
will  start  to  propagate  sooner  than  in  the  absence  of  transformation  because 
AK  is  initially  positive  due  to  AKg.  As  the  crack  grows,  AK  diminishes 
and  eventually  becomes  negative  at  about  Aa/w  =0.1.  This  means  that  the 
applied  loads  must  be  increased  to  sustain  crack  growth  and  in  a  stiff  loading 
system  the  crack  will  propagate  stably  under  rising  load.  This  will  continue 
until  Aa/w  is  about  0.7  and  the  SIF  change  is  then  given  by 


(1-v)  AK 
E  eT  vf  ^ 


0.35. 


(23) 


This  expression  represents  the  toughness  enhancement  and  is  larger  than  that 
given  by  equation  (18).  After  this  amount  of  crack  growth,  AK  increases  and 
so  if  the  loads  are  kept  constant  or  Increased,  the  applied  K  will  exceed 
Kjq.  Thereafter  the  propagation  will  become  unstable.  Thus  the  maximum 
magnitude  of  AK  represented  by  equation  (23)  is  equivalent  to  the  asymptotic 
value  of  AK  due  to  dilatation  alone  observed  by  McMeeking  and  Evans  [19]. 

It  is  interesting  to  note  that  the  asymptotic  value  of  AK  due  to  dilatation 
plus  shear  is  not  the  relevant  quantity.  The  useful  toughness  enhancement  is 
due  to  the  minimum  in  the  curve  for  AK  at  Aa/w  =  0.7.  Thus  it  may  not  be 


been  extended  to  the  cases  when  X  *  0  . 

In  contrast  with  AKq,  AKg  Is  not  zero  and  it  does  depend  on  the  value 
of  X.  In  fact,  AKg  Increases  as  X  Increases.  Equation  (16)  indicates 
that  the  shearing  effect  results  in  fracture  toughness  reduction  because  AKg 
is  positive.  This  reduction  for  X  -  1  is  comparable  to  the  fracture 
toughness  enhancement  in  the  case  of  a  purely  dllatant  transformation  given  by 
equation  (18)  when  the  crack  has  advanced  Aa  *  5w  .  The  AKg  results  in 
table  1  can  almost  be  reproduced  by  the  equation  (17)  for  the  respective 
values  of  X.  This  means  that  the  transformed  material  and  the  shearing 
contribution  do  not  affect  the  features  of  the  transformation  zone  as  it  is 
found  from  the  unperturbed  elastic  crack  tip  field.  Therefore,  the  fact  that 
AKg  is  positive  is  a  consequence  of  the  nature  of  the  equation  (16).  An 
attempt  to  justify  this  may  be  made  by  regarding  the  range  of  the  positive 
contribution  to  the  Integral  of  equation  (16).  Since  for  the  stationary  crack 

3Q 

Q  ■  ■£—  ,  the  integral  sign  depends  on  the  sign  of  the  integrand  sin(20). 

Hence,  it  is  positive  when  9  <  ^  and  negative  when  9  >  y  .  Therefore,  it 
can  be  said  that  the  positive  contribution  comes  from  a  large  sector  and  can 
very  likely  override  the  negative  contribution. 

The  consequence  of  the  computed  AKg  values  would  be  that  crack 
propagation  may  take  place  sooner  than  when  there  is  no  transformation.  This 
is  because  AKg  is  positive  for  the  stationary  crack  and  KA  can  become 
greater  than  Kj£  earlier  than  otherwise.  However,  this  growth  is  likely  to 
be  stable,  a  point  which  is  elucidated  in  what  follows. 

Propagating  crack 

In  order  for  the  crack  to  advance  quasi-statically  it  is  required  that 

A  A 

the  applied  loads  KA  be  such  that  KA  ■  Kj£  where  Kj£  is  given  by 


equation  (16).  Therefore,  by  knowing  AK  ,  we  can  produce  an  "R  curve"  [ 1 5 J 


-C.  vVi  -  V- 
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zones  are  developed  Inside  a  circle  centered  at  the  crack  tip  whose  area  Is 
less  than  0.15Z  of  the  mesh.  Indeed,  these  nuabers  fall  within  the 
corresponding  length  range,  set  for  saall  scale  yielding  condition  In  the 
elastic-plastic  fracture  analysis  by  Rice  and  Tracey  [23].  In  addition,  the 
fact  that  the  displacements  of  the  elements  near  the  perimeter  of  the  domain 
(fig.  5)  are  the  unperturbed  elastic  crack  tip  field  displacements,  ensures 
the  small  scale  condition  too.  In  figure  9,  it  can  be  seen  that  the  zone 
taken  numerically  for  X  *  0  is  almost  identical  with  that  used  by  McMeeking 
and  Evans  [19]  as  predicted  by  Budlansky  et  al.  [20].  This  means  that  the 
transformed  material  does  not  affect  the  linear  elastic  crack  tip  field  given 
by  equation  (5).  The  same  argument  applies  to  the  case  when  X  *  1  (fig. 

11).  However,  the  shearing  contribution  affects  the  zone  shape  especially  at 
small  angles  of  9.  It  is  also  worth  mentioning  that  the  zone  height  w  is 
the  same  for  all  values  of  X  as  shown  in  table  1. 

Consider  now  the  SIF  change  AKq.  It  can  be  deduced  from  table  1  that 
AKp  -  0  if  we  take  into  account  the  numerical  error  involved  in  the 
calculations.  Therefore,  for  X  >  0  no  AKD  is  observed  even  though  the 
zone  boundaries  differ  from  that  with  X  ■  0.  This  independence  of  AKD  from 

X  can  be  justified  by  observing  equation  (13).  The  sign  of  the  integral 

36  w 

depends  on  the  sign  of  the  cosOj— ).  The  integral  is  positive  for  6  <  y  and 

negative  for  0  >  y  .  Thus,  the  positive  contribution  to  the  integral  value 

comes  from  a  much  smaller  section  than  the  sector  of  the  negative 

contribution.  It  is  this  wider  ranging  source  of  the  negative  effect  which 

dominates  the  integral  value.  Therefore,  the  larger  range  of  r  for  small  6 

when  X  deviates  from  zero  cannot  be  cause  of  a  significant  disturbance  of 

the  final  result.  This  remark  has  already  been  made  before  [19],  but  only 

within  the  frame  of  a  purely  dllatant  transformation.  The  argument  now  has 


The  crack  has  been  assumed  to  propagate  quasl-statlcally  under  the 
conditions  mentioned  in  section  4.  The  SIF  change  calculations  AKD  and 
AKg  were  carried  out  numerically  for  a  given  crack  advance  Aa.  The  formulae 
(13)  and  (14)  were  used  again*  The  results  are  shown  in  tables  2  through  4. 

In  figures  12  and  13  the  changes  AKp  and  AKg  are  plotted  against  the  crack 
advance  Aa  for  X  *  0  and  X  -  1  respectively*  In  table  5  the  symptotic, 
l*e.  the  maximum  SIF  change  AK  *  AKp  +  A  Kg  is  shown. 

8.  DISCUSSION 

In  this  discussion  of  the  results  quoted  in  the  previous  section,  we 

shall  focus  mainly  on  how  the  shear  component  of  the  transformation  influences 

the  fracture  toughness  behavior  of  the  material.  It  should  be  borne  in  mind 

T  T 

that  the  shearing  contribution  is  characterized  by  the  parameter  X  •  y  /e  • 
Large  values  of  X  denote  large  transformation  shear  strains  for  a  given 
volume  dilatation. 


Stationary  crack 


It  has  been  mentioned  before  that  the  dllatant  part  of  the  transformation 


does  not  affect  the  material  fracture  toughness,  l.e.  AKp  -  0.  Furthermore, 
Lambropoulos  [21]  based  on  a  transformation  zone  derived  from  the  unperturbed 


linear  elastic  crack  tip  field,  concluded  that  AKg  »  0  too.  Our  analysis' 
predictions  are  in  accord  with  the  above  results  only  as  far  as  AKp  is 
concerned.  Table  1  shows  that  AKg  is  not  zero.  However,  we  have  used  a 
different  constitutive  law  from  Lambropoulos. 


Before  discussing  the  nature  of  AKp  and  AKg  we  should  emphasize  that 


the  small  scale  transformation  condition  is  satisfied.  This  Is  because  the 


Stationary  crack 

The  transformation  gone  for  A  m  0  la  shown  In  fig.  8.  There  Is  no 
shearing  effect  because  ■  0.  The  asterisks  Indicate  the  element 
Integration  stations  where  the  condition  for  transformation  was  met.  The 
transformed  zone  covers  an  area  only  0.1Z  of  the  mesh.  The  nodal  displace¬ 
ments  of  the  elements  near  the  external  boundary  of  the  domain  (fig.  S)  are 
Identical  within  3  significant  figures  to  the  linear  elastic  ones  when 
transformation  does  not  take  place.  The  periphery  of  the  transformation  zone 
Is  shown  in  fig.  9.  The  curve  denoted  by  the  letter  D  is  the  periphery  of  the 

transformation  zone  used  by  McMeeking  and  Evans  [19].  This  last  zone  is  given 
8  2  6 

by  r  * - w  cos  (•*■)  where  the  height  w  has  been  taken  to  agree  with  the 

3/3 

height  of  the  finite  element  zone. 

The  transformation  zone  for  X  ■  1  is  shown  in  fig.  10.  The  boundary  of 
the  zone  is  shown  in  fig.  11  along  with  the  zone  used  by  McMeeking  and  Evans 
as  described  before.  The  zones  have  the  same  height  as  the  zone  for  A.  *  0. 
However,  its  leading  front  appears  to  be  elongated  along  the  x  axis  as  y  >  0 
compared  to  the  other  zone.  As  a  result,  an  inflection  point  is  observed. 

The  zone  boundary  to  the  left  of  this  point  tends  to  be  concave  whereas  the 
boundary  to  the  right  tends  to  be  convex.  It  is  worth  noting  that  this 
tendency  for  elongation  increases  with  increasing  values  of  X.  This 
observation  has  been  confirmed  by  the  appearance  of  the  zone  boundaries  that 
we  generated  for  intermediate  values  of  X. 

The  SIF  change  calculations  were  carried  out  numerically  by  using  the 


formulae  (13)  and  (14)  for  the  AKD  and  AKg  changes  respectively.  These 
changes  are  shown  in  Table  1. 


where  [C]  Is  Che  linear  elastic  stress-strain  matrix  for  plane  strain,  for  the 
composite,  (T)  is  the  two-dimensional  traction  vector,  {u^}  is  the  nodal 
displacements  and  {E^,}  is  the  element  transformation  strain  vector  given  by 

eT  -  yt  sln2Q 
3(l~2v)  6  Y  8lnZU 

(et1  -  eT  +  yT  8in2Q  (22) 

yT  cos2Q 

It  is  worth  noting  that  the  stiffness  matrix  [K]  is  independent  of  the 
transformation  since  it  involves  only  the  elasticity  tensor  [C]  and  the 
Interpolation  matrix  [B].  The  same  is  true  for  the  vector  {F^}  .  However, 
the  quantity  {F^}  depends  on  the  transformation  through  {E^}  . 

An  Iterative  method  was  used  to  solve  equations  (21)  as  stated  in  section 
5,  All  calculations  were  carried  out  for  Poisson's  ratio  v  -  0.25.  Solutions 
were  obtained  for  X  equal  to  0,  0.25,  0.5,  0.75  and  1.  Extremely  rapid 
convergence  was  achieved  with  3  iterations  when  the  computed  displacements 
were  found  to  converge  to  within  3  significant  figures. 

7.  RESULTS 

In  this  section  we  shall  discuss  the  results  of  the  finite  element 
analysis  for  the  stationary  crack  problem.  We  shall  also  give  the  SIF  changes 
for  both  stationary  crack  and  quasi-statically  advancing  crack  deduced  through 
a  model  for  that  case.  The  shape  and  features  of  the  transformation  zone  for 
the  propagating  crack  are  determined  as  discussed  in  section  4  starting  from 
the  stationary  crack  zone.  However,  the  stationary  crack  zone  is  now  the  one 
deduced  via  the  finite  element  calculations. 
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[  EkJt  Cijkl  6eij  ^  "  £  Ti  6ui  d*  +  [  ekA  Cljki  6eiJ  ^  (20) 

Equation  (20)  illustrates  the  nonlinearity  of  the  problem.  That  is,  the 
second  integral  on  the  right  hand  side  can  only  be  determined  after  the 
solution  is  found*  Thus,  the  solution  procedure  will  Involve  iteration. 
Assuming  a  solution  we  determine  the  transformation  zone.  Next,  we  solve  the 
resulting  linear  problem  arising  from  equation  (20)  with  the  second  Integral 
on  the  right  hand  side  evaluated  and  compare  the  solution  with  the  assumed 
one.  The  process  continues  until  convergence  is  achieved.  The  linear  elastic 
solution  can  be  used  to  initiate  the  process  with  eA  ■  0  i.e.  with  no 
transformed  region.  This  method  is  addressed  in  detail  in  the  next  section. 

6.  FINITE  ELEMENT  CALCULATIONS 

The  finite  element  method  was  used  to  solve  the  boundary  value  problem  as 
stated  in  the  previous  section.  The  domain  (fig.  5)  was  discretized  into  554 
4-noded  quadrilateral  isoparametric  elements  with  4  integration  stations.  In 
figures  6  and  7  the  224  element  near  tip  mesh  and  the  320  element  far  field 
mesh  are  shown  respectively.  The  second  mesh  surrounds  the  first  one.  The 
fine  near  tip  discretization  was  used  in  order  to  ensure  that  the 
transformation  zone  Is  determined  fairly  reliably.  By  using  the  standard 
element  displacement  and  strain  interpolation  matrices  [A]  and  [B] 
respectively,  we  rewrite  the  governing  equation  (20)  as  follows 

IK]  {uj,}  -  {Fb)  +  {Ft>  (21a) 

with  [Kj  -  /  [B]T  IC]  [Bj  dA  (stiffness  matrix)  (21b) 

A 

{F.}  -  /  [AJT  (T)  ds  (applied  load  vector)  (21c) 

0  A 

(Ft>  ■  /  [B]^  {E,j}  ds  (transformation  load  vector)  (21d) 
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transformation  zone  appropriately*  This  is  accomplished  by  taking  into 
account  the  perturbation  of  the  linear  elastic  crack  field  induced  by  the 
transformation*  This  subject  is  addressed  in  the  next  section. 

5.  THE  BOUNDARY  VALUE  PROBLEM  FOR  THE  STATIONARY  CRACK 

We  are  concerned  with  the  stresses  and  deformations  near  the  tip  of  a 
long  crack  in  plain  strain.  The  body  is  loaded  so  that  only  mode  I  (tensile 
opening  mode)  stresses  arise  at  the  tip.  We  are  concerned  with  "small  scale 
transformation"  where  the  transformation  zone  is  confined  to  a  region  very 
close  to  the  crack  tip.  This  can  be  achieved  by  imposing  tractions  on  a 
circular  boundary  far  from  the  tip  in  agreement  with  the  standard  singular 
linear  elastic  solution  given  by  equation  (5).  Symmetry  permits  the  analysis 
to  be  confined  to  a  semi-circular  domain  as  shown  in  fig.  5.  The  crack 
surface  is  traction  free  and  the  symmetry  line  is  free  of  shear  traction  and 
displacements  normal  to  the  line. 

The  governing  equations  in  the  plane  can  be  stated  by  the  principle  of 
virtual  work  in  the  absence  of  body  forces. 

I  o  6e±j  dA  -  /  T±  6Ul  ds  (19) 

where  A  is  the  area  and  S  is  the  perimeter  of  the  domain,  T^  the 
tractions  on  S  and  u^  the  displacements  in  A.  The  symbol  6  indicates 
an  arbitrary  virtual  variation  of  the  quantity  it  precedes.  In  addition  the 
stresses  are  related  to  the  strains  by  the  constitutive  law  given  by  equation 
(2).  As  a  result,  the  governing  equation  (19)  can  be  rewritten  as  follows 


function  of  Aa  ,  the  amount  of  crack  propagation*  The  slope  versus  Aa 
decreases  continuously  and  the  curve  appears  to  approach  an  asymptote  for 
large  values  of  Aa*  At  that  stage  [19,20] 

(1-v)  AK_ 

- - - 3  -  -  0.22  (18) 

E  e  Vj  rw 

To  maintain  the  critical  propagation  value  of  K  ■  KIC  for  the  material  at 

the  crack  tip  the  loads  would  have  to  be  increased  with  Aa.  This  will  mean 

an  apparent  higher  value  of  K  «  computed  from  the  applied  loads.  The 

material  will  appear  to  be  toughened  as  a  result. 

Consider  now  the  shearing  contribution  to  AK.  According  to  our 

hypothesis,  particles  enter  the  transformation  zone  ahead  of  the  crack  and  the 

direction  of  shearing  transformation  will  be  determined  by  the  state  of  stress 

there.  We  assume  still  that  the  stress  field  is  unperturbed  by  the  presence 

30 

of  a  zone  and  so  Q  equals  to  current  value  of  measured  from  the  crack 

tip.  As  the  crack  propagates  by,  this  value  of  Q  will  remain  unchanged  for 
that  particle.  Thus  all  particles  on  a  line  parallel  to  the  crack  will  have 
the  same  value  of  Q  ,  except  those  in  the  zone  created  before  the  crack 
propagated  (fig*  3).  The  resulting  contribution  AK§  from  such  a  wake  zone 
is  negative.  The  maximum  of  the  absolute  value  of  AK<*  for  all  X.  is 
achieved  when  the  crack  propagates  0.5w  whereas  the  maximum  of  the  absolute 
value  of  AKq  as  found  from  equation  (18),  is  achieved  at  Aa  **  5w.  In 
figure  4  the  ratio  | |  max ^ I  Climax  *8  8hown  plotted  against  the  ratio  X. 

As  it  is  seen,  the  amount  of  the  shearing  contribution  to  the  toughening  can 
become  the  same  as  that  of  the  dilatant  contribution  when  X  approaches  1. 

In  conclusion,  the  importance  of  the  shearing  contribution  raises  the  need  for 
a  more  rigorous  estimation.  Therefore,  it  is  required  that  we  compute  the 


Results  for  stationary  crack 

In  this  case,  the  transformed  zone  has  a  perimeter  which  Is  the  locus  of 
points  of  equal  hydrostatic  'tress*  McMeeklng  and  Evans  [19]  and  Budlansky  et 
al.  [20]  have  shown  that  AKD  ■  0  for  this  case. 

As  far  as  the  shearing  contribution  Is  concerned,  Q  Is  taken  to  be  the 
angle  to  the  principal  shearing  direction  when  the  material  transformed,  l.e. 
at  the  edge  of  the  zone  then  prevailing.  The  stress  state  Is  taken  to  be  the 

30 

unperturbed  linear  crack  tip  field  and  thus  Q  *  -y-  •  As  a  result 

(1-v)  AK„ 

- = - ~  ■  0. 15X  (17) 

E  e  vf  /w 

where  w  is  the  height  of  the  transformation  zone  as  shown  in  fig.  2,  X  is 
the  ratio  yT/eT  and  v  »  0.3.  The  transformation  strain  has  been  rewritten 
as  Vf  Yp  where  Yp  is  the  unconstrained  transformation  strain  after 
twinning  for  an  individual  particle. 

Equations  (16)  and  (17)  indicate  that  the  shearing  effect  causes  fracture 
toughness  reduction.  This  reduction  becomes  larger  as  X  increases  i.e.  the 
shear  strain  yT  increases  for  a  given  e^.  As  a  consequence,  crack 
propagation  may  occur  more  readily  than  in  the  absence  of  transformation. 


Results  for  propagating  crack 

If  the  crack  propagates  quasi-statically  and  there  is  no  reverse 
transformation,  a  wake  of  transformed  material  is  left  along  the  newly  created 
crack  surfaces.  This  wake  will  be  parallel  to  the  crack  surfaces  if  we  assume 
that  the  tip  propagates  at  a  constant  value  of  R  (fig.  3).  The  dilatant 
contribution  to  AK  due  to  such  a  transformed  zone  has  been  studied  by 
McMeeklng  and  Evans  [19]  and  has  been  found  to  be  a  monotonically  decreasing 
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McMeeking  and  Evans  [19]  have  shown  that 


“d— Jr -  fr 

6  /2n  (1-v)  A^ 


3 
"  7 


,30. 

cos(-j— )  dA 


(13) 


where  E  is  the  Young's  modulus*  From  the  form  of  the  weight  function  h^ 

T 

and  e  it  can  be  deduced  that 

T  -  — 

AKC - 3E  ?■  Jr2  cos(2Q  -  £■)  sinG  dA  (14) 

8  /ST  (l-vZ)  A^, 

30 

The  fact  that  cos(-^— )  is  an  even  function  of  0  allows  the  integral  in 
equation  (13)  to  be  computed  in  the  region  0  <  0  <  u  (figures  2,3).  The 
same  argument  applies  to  equation  (14)  because  yT  sin0  is  an  even  function 
of  9.  That  is  because  both  yT  and  sin0  are  odd  functions  of  0. 

The  SIF  change  AK  is  related  to  the  applied  field  KA  as  follows 

1^.  -  KA  +  AK  (15) 

where  Kt  is  the  S.I.F.  that  actually  describes  the  crack  tip  field  when  the 
transformation  has  taken  place.  If  AK  is  negative  equation  (15)  Implies 
that  the  stress  intensity  at  the  crack  tip  is  lower  than  that  associated  with 
the  applied  loads.  In  particular,  since  the  crack  propagates  when  Kj 
becomes  equal  to  the  fracture  toughness  KIC  of  the  pretransformed  material 
[14]  the  apparent  fracture  toughness  Kjq  is  given  by 

-  *10  -  “  (16) 

If  AK  is  negative  the  material  appears  to  be  tougher  than  the  pre transformed 


one. 
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TABLE  CAPTIONS 


Table  1.  Stress  Intensity  factor  changes  for  the  stationary  crack  calculated 
by  the  finite  element  method. 

Table  2.  Stress  Intensity  factor  changes  for  the  propagating  crack  when 
X  -  *  0  deduced  from  the  finite  element  calculations. 


Table  3.  Stres^  intensity  factor  changes  for  the  propagating  crack  when 
\  m  y  /£  m  o.5  deduced  from  the  finite  element  calculations. 


Table  4.  Stresjj  intensity  factor  changes  for  the  propagating  crack  when 
X  *  y  /e  -I  deduced  from  the  finite  element  calculations. 

Table  5.  Maximum  fracture  toughness  changes  achieved  after  the  crack  has 
advanced  quasl-statically. 


Table  1 


Stress  Intensity  Factor  Changes  for  the  Stationary  Crack 

(1-v)  AKjj  (l-v)  AKg 


w 

(zone  height) 

E 

T  j— 
e  rw 

T 

E  e1  Vj  i 

0.00 

2.14 

-0.0085 

0.0000 

0.25 

2.14 

-0.0066 

0.0381 

0.50 

2.14 

—0.0066 

0.0762 

0.75 

2.14 

-0.0049 

0.1145 

1.00 

2.14 

-0.0049 

0.1527 

akd  - 

Dilatational 

SIF 

change 

AKS  = 

Shearing 

SIF 

change 
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Table  2 


Stress  Intensity  Factor  Changes  for  the  Advancing  Crack  When  X  ■  0 


A a 

(1-v)  AKp 

(1-v)  AKg 

(1-v)  AK 

w 

T  — 

E  e  Vj  /w 

E  eT  vf  ^ 

E  e  /w 

0.0 

-0.0085 

0 

-0.0085 

-0.0780 

0 

-0.0780 

0.5 

-0.1270 

0 

-0.1270 

0.8 

-0.1550 

0 

-0.1550 

1.1 

-0.1723 

0 

-0.1723 

2.0 

-0.1961 

0 

-0.1961 

3.0 

-0.2055 

0 

-0.2055 

4.0 

-0.2097 

0 

-0.2097 

5.0 

-0.2119 

0 

-0.2119 

7.0 

-0.2142 

0 

-0.2142 

10.0 

-0.2157 

0 

-0.2157 

15.0 

-0.2167 

0 

-0.2167 

-0.2171 

0 

-0.2171 

30.0 

-0.2174 

0 

-0.2174 

Aa 

■  Dimensionless 

crack 

.  advance 

w 

akd 

*  Dilatational 

SIF 

change 

AKg 

*  Shearing 

SIF 

change 

Total 


SIF  change 


Table  3 


Stress 

Intensity  Factor  Changes  for 

the  Advancing 

Crack  When  X.  -  0.5 

Aa 

(1-v)  AKjj 

(1-v)  AKg 

(1-v)  AK 

w 

E  e^  /w 

E  e^  Vj  /w 

T 

E  e  vf  ✓ 

-0.0066 

0.0762 

0.0695 

-0.0961 

-0.0954 

-0.1915 

0.6 

-0.1362 

-0.1057 

-0.2419 

1.0 

-0.1656 

-0.0860 

-0.2516 

2.0 

-0.1943 

-0.0376 

-0.2319 

3.0 

-0.2037 

-0.0312 

-0.2348 

4.0 

-0.2079 

-0.0298 

-0.2377 

5.0 

-0.2101 

-0.0287 

-0.2388 

10.0 

-0.2139 

-0.0268 

-0.2407 

15.0 

-0.2149 

-0.0258 

-0.2406 

20.0 

-0.2153 

-0.0247 

-0.2400 

30.0 

-0.2156 

-0.0230 

-0.2386 

Aa 

w 

■  Dimensionless  crack 

advance 

akd 

■  Dilatational  SIF 

change 

AKS 

■  Shearing  SIF 

change 

SIF  change 


Total 


(1-v)  AKp 


(1-v)  AKg 


(1-v)  AK 


E  e  vf  fit 

E  e  /w 

E  e1  vf  /w 

-0.0049 

0.1527 

0.1479 

-0.0944 

-0.1904 

-0.2848 

-0.1344 

-0.2110 

-0.3454 

-0.1638 

-0.1718 

-0.3356 

-0.1925 

-0.0744 

-0.2669 

-0.2019 

-0.0620 

-0.2638 

-0.2061 

-0.0592 

-0.2653 

-0.2084 

-0.0570 

-0.2653 

-0.2121 

-0.0532 

-0.2653 

-0.2131 

-0.0512 

-0.2643 

-0.2135 

-0.0492 

-0.2626 

-0.2138 

-0.0455 

-0.2594 

Dimensionless  crack  advance 

Dilatatlonal  SIF  change 

Shearing  SIF  change 

Total  SIF  change 

Table  5 


SIF  change 
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Pig.  1 

Fig.  2. 

Fig.  3. 

Fig.  4. 

Fig.  5. 

Fig.  6. 
Fig.  7. 

Fig.  8. 

Fig.  9. 
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Figure  Captions 


a.  Direction  Q  of  the  maximum  shear  stress  i 

b.  Shape  change  of  an  unconstrained  transforming  element.  The 
shearing  direction  is  along  the  maximum  shear  stress  direction. 


Transformation  zone  shape  based  on  critical  hydrostatic  stress 
transformation  and  the  unperturbed  elastic  solution  for  the 
stationary  crack.  The  maximum  height  w  of  the  zone  appears  at  $ 
°  60°.  The  shearing  direction  Q  at  transformation  is  also 
Indicated  at  an  arbitrary  point  (r,0). 


Transformation  zone  shape  for  the  quasi-statically  advancing  crack. 
The  flank  of  the  zone  is  parallel  to  the  crack  surface  and  tangent 
to  the  leading  fro  at  of  the  zone  at  8  ■  60®.  The  shearing 
direction  of  the  newly  transforming  elements  is  determined  by  the 
angle  Q  at  the  leading  front  of  the  zone.  The  initially 
transformed  elements  retain  their  own  shearing  directions. 


Comparison  of  the  maximum  shearing  contribution  to  the  toughness 
enhancement  with  the  maximum  dilatant  contribution.  Those 
contributions  are  derived  from  a  zone  shape  shown  in  fig.  3  when 
the  crack  propagated.  The  zone  for  the  stationary  crack  is  shown 
in  fig.  1.  The  AKS  is  achieved  after  Aa  »  0.5w  whereas  the 
AK0  is  achieved  after  Aa  »  5w. 


Domain  and  boundary  conditions  for  the  boundary  value  problem  for 
the  stationary  crack. 


Near  tip  finite  element  mesh. 


Far  field  finite  element  mesh  which  surrounds  the  mesh  shown  in 
fig.  6. 


Transformation  zone  for  the  stationary  crack  in  the  case  of  purely 
dilatant  transformation  (y  -  0)  derived  by  solving  the  boundary 
value  problem  by  the  finite  element  method. 


Comparison  of  the  zone  boundary  of  fig.  8  with  the  zone  boundary 
derived  by  McMeeking  and  Evans  (19]  for  the  purely  dilatational 
transformation  and  marked  D.  (fig.  2). 
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Fig.  10.  Transformation  zone  for  the  stationary  crack  when  the  unconstrained 
shear  transformation  strain  is  equal  to  the  unconstrained  volume 
dilatation.  The  zone  has  been  derived  by  solving  the  boundary 
value  problem  by  the  finite  element  method. 


Fig.  11.  Comparison  of  the  zone  boundary  of  fig.  10  with  the  zone  boundary 
derived  by  McMeeking  and  Evans  [19]  for  the  purely  dilatational 
transformation  (fig.  2)  (D). 

Fig.  12.  The  R-curve  predicted  from  the  stress  Intensity  factor  analysis  in 
the  case  of  purely  dilatational  transformation. 

Fig.  13.  The  R-curve  predicted  from  the  stress  intensity  factor  analysis  in 
the  case  when  the  unconstrained  shear  transformation  strain  is 
equal  to  the  unconstrained  volume  dilatation. 
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